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Abstract 

Principal Hugoniot and K-shell X-ray absorption spectra of warm dense KC1 are calculated using 
the first-principles molecular dynamics method. Evolution of electronic structures as well as the 
influence of the approximate description of ionization on pressure (caused by the underestimation 
of the energy gap between conduction bands and valence bands) in the first-principles method are 
illustrated by the calculation. Pressure ionization and thermal smearing are shown as the major 
factors to prevent the deviation of pressure from global accumulation along the Hugoniot. In 
addition, cancellation between electronic kinetic pressure and virial pressure further reduces the 
deviation. The calculation of X-ray absorption spectra shows that the band gap of KC1 persists 
after the pressure ionization of the 3 p electrons of Cl and K taking place at lower energy, which 
provides a detailed understanding to the evolution of electronic structures of warm dense matter. 
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I. INTRODUCTION 


Warm dense matter (WDM) is of particular interest to astrophysics, geophysics, and in¬ 
ertial confinement fusion (ICF).[l-3] Recent developments of high-power laser facilities^, 5| 
and corresponding diagnosing techniques [6], 7] greatly facilitate its investigation, making 
WDM accessible to systematical laboratory studies. Laser-driven shock compression is a 
typica 


IIU 


eV 


way to generate warm dense states in large-scale laser facilities, such as Shen-Guang 
9fl. Among materials used in shock compression experiments, KC1 is a good choice 


for X-ray absorption measurements) 


considering its large intrinsic band gap of 8.69 


111 ]. which efficiently reduces the preheating of diagnosing X-ray in the early stages of 


the experiments 


10], As a typical ionic crystalline material (under ambient conditions) 


composed of heterogeneous species, KC1 is also beneficial for illustrating the evolution of 
electronic structures (e.g., pressure ionization, energy level shifts, and metalization) of het¬ 
erogeneous ionic materials under warm dense conditions. Those properties can now be 
readily detected by their X-ray absorption spectra (XAS). In particular, the variation of the 
intrinsic band gap under compression has a significant influence on the degree of ionization, 
which is a critical parameter in several commonly used models of equation of state (EOS), 
e.g., the quotidian equation of state (QEOS) model 12], 

From the methodological point of view, KC1 is also an illustrating example to display the 
influence of the approximate description of ionization, caused by the underestimation of the 
energy gap between conduction bands and valence bands, in the first-principles molecular 
dynamics (FPMD) method, now an influential theoretical tool in the study of WDM [r]]. 
The FPMD method, usually implemented with local density approximation (LDA) or the 
Perdew-Burke-Ernzerhof (PBE) version of generalized gradient approximation (GGAlfpfj] 
to the exchange-correlation interaction, has been known to underestimate the energy gap 
by 40% - 60% for nonmetals[ljj], which causes a serious overestimation in the ionization, 
and was thus postulated to bring about sizable deviations in the calculation of thermal 
properties of WDM. However, this conjecture is not well supported by recent studies. The 
deviation turns out to be small. FPMD calculations on a large variety of materials, including 
hydrogen, 16ql8] helium, 19j aluminum, 20,121] and iron 22], display a good agreement with 
measured Hugoniot data. These calculations cover a large range of temperature and pressure, 
at which electrons in the outermost shell and inner shell can be ionized. How to understand 
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this small effect remains an open problem. The answer to it, however, is important to the 
improvement of the FPMD method. The ionization of KC1 takes place under a condition 
relatively easy to realize through current WDM experimental facilities, which makes it a 
favorable example to examine this effect. 

Using the FPMD method, we investigate the variation of electronic structures and the 
influence of the approximate description of ionization to the calculation of pressure in warm 
dense KC1. The variation of electronic structures is illustrated through the calculation of 
XAS. Our results show that the energy gap, originally between the occupied valence band 
and unoccupied conduction band under ambient conditions, persists after the occurrence of 
pressure ionization of the 3 p electrons of Cl and K taking place at a lower energy. Pressure 
ionization and thermal smearing are shown as the major factors to prevent the deviation 
of pressure caused by the approximate description of ionization from global accumulation 
along the Hugoniot. In addition, cancellation between electronic kinetic pressure and virial 
pressure further reduces the deviation, leading to a small influence of the approximate de¬ 
scription of ionization in the FPMD method, in line with the trend illustrated in previous 


calculations 


ld-|22|. 


The rest of the article is organized as follows. In Sec. II methods and numerical details 
used in the calculation of principal Hugoniots and K-shell XAS are summarized. In Sec. 
Ill, the influence of the approximate description of ionization on the calculation of pressure 
is discussed through the calculation of principle Hugoniot. In Sec. IV, the variation of 
electronic structures of KC1 in warm dense states is displayed together with the calculation 
of XAS. We conclude the article with a short summary in Sec. V. 


II. COMPUTATIONAL METHODS AND NUMERICAL DETAILS 


A. Calculation of Principal Hugoniots 


Principal Hugoniot of shock-compressed KC1 is calculated through FPMD simulations 
consisting of 54 atoms, i.e., 27 formula units of KC1. Finite-temperature density functional 
theory (FTDFT)Q is employed in the simulation to describe the statistics of electrons. 


The calculation is performed using the Quantum-Espresso package. 24] The PBE version of 
GGaQ is used to account for the exchange-correlation interaction. Projector augmented 
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wave (PAW) pseudopotentials[25] generated by the ATOMPAW program 26] are employed to 
model the electron-ion interactions. The plane-wave cutoff is set to be 50 Ry, which controls 
the deviation within 0.01 eV in the calculation of total energy. Only the T point is used to 
sample the Brillouin zone in the FPMD simulations. 

For a given density, a series of simulations are carried out at different ionic temperatures, 
which are controlled by the Nose-Hoover thermostat 27]. The electronic states are populated 
according to the Fermi-Dirac distribution and thermodynamical equilibrium is maintained 
by setting the electronic temperature (T e ) the same as the ionic temperature (7)). In some 
cases, a T e different from the value of T) is used to illustrate the influence of the approximate 
description of ionization. Typical dynamic simulations last for 3 ~ 5 ps with an appropriate 
time step 0.4 ~ 1.0 fs, varying with temperature and density. Data presented in the article 
are collected from the last 2000 steps. 

The principal Hugoniot points are determined by solving the Hugoniot equation 28] 


£i-£ 0 = -(Pi + Po)(Ro-Ri), 


( 1 ) 


where E represents the average internal energy, V is the average volume, and P is the 
pressure. Subscripts 0 and 1 refer to initial and shocked states, respectively. Kinetic pressure 
of ions, equal to k B T/V, is added manually in the calculation to accurately account for the 
total pressure. The reference state (E 0 , Vo) is determined from a face-centered cubic (fee) 
structure under ambient conditions. Pq is neglected in the calculation as it is several orders 
smaller than P\. An auxiliary cubic polynomial interpolating function is employed in solving 
Eq. ([T]) to reduce the numerical error less than 1%. 

Note that KC1 undergoes a polymorphic phase transition from a fee lattice (Bl, NaCl 
structure) to a body-centered cubic (bcc) lattice (B2, CsCl structure) during the shock 


compression 


29 


32], T 


re transition pressure varies from 1.97 to 2.5 GPa, depending on 


experimental conditions 


29 


32], The calculation of Hugoniot starts from a B2 structure at 


the pressure around 3.5 GPa, and phase transition points along the Hugoniot are avoided. 


B. Calculation of XAS 


The XAS of warm dense KC1 is calculated via first-principles methods based on FPMD, 


similar to those used in previous studies. 


33 


34] Theoretically, X-ray absorption is charac- 
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terized by its cross-section a(u), which is calculated in a perturbative way as [35 1 

<J ( UJ ) 


47 r 2 ahu 


= EK/I 

f 


e • r 




( 2 ) 


where, letters / and % represent final and initial states respectively, e is a unit vector de¬ 
noting the polarization of incident light, a is the fine structure constant, and E(Ef) is 
the Fermi-Dirac distribution of final states at finite temperature. Calculated using Eq. ([2]), 
K-shell XAS of shock-compressed A1 have displayed a good agreement with experimental 

lermodynamical properties 


measurements 


Ml, 


36]. Information on electronic structures and_ 


of WDM can be further derived through these calculations. [341, 


36] 


The XAS calculation is performed using the XSPECTRA package 35, |37| supplied with 


the Quantum-Espresso distribution 24], A minor modification is made to account for the 
Fermi-Dirac distribution of electrons at finite temperature. Wave functions of the final 
states used in the XAS calculation are prepared by a density functional theory (DFT) 
calculation on the atomic configurations generated in the previous Hugoniot calculation. 
The plane-wave energy cutoff is 50 Ry, and the hole in the Is state of the excited Cl 
atom is manually included in the core of the pseudopotential with the help of the GIPAW 


pseudopotential 


38] to approximately account for electron-hole interactions. Other atoms 


m 


the system are presented using the Troullicr-Martins pseudopotentials 39]. The exchange- 
correlation functional is kept the same as that used in the Hugoniot calculation. With the 
wave functions prepared, cr(uj) is then calculated through Eq. ([2]). For each thermodynamical 
condition, a total of 8 XAS are calculated on selected snapshots along the trajectory of ions, 
separated by 200 time steps. The final results of XAS are averages of the spectra of each 
snapshot. 

The position of the K-absorption edge, defined as E edge - Ei s , is calculated in two steps. 
Firstly, the edge energy with respect to the chemical potential /q E edge - /q is determined 
directly from the XAS as the intersection of the abscissa to the slope of the edge. Secondly, 
the energy of the Is orbital of Cl referring to /q i.e., E\ s - /q is calculated by an separate 
FPMD simulation on a system consisting of 16 atoms. 8 atomic configurations separated 
by 200 time steps are generated by the FPMD simulation. Then, for each configuration, a 
DFT calculation with pseudopotentials explicitly including all electrons (generated by the 
AT0MPAW pseudopotential package [40] with a core cutoff radius of 0.5 Bohr and a plane- 
wave cutoff energy of 400 Ry) is performed to get Ei s - q. The difference between these two 
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FIG. 1. (Color online) Shocked principal Hugoniots of KC1 in (a) P -p diagram and (b) T -p diagram. 

energies is the K-edge position E edge - E ls required. 

Note that in the second step, the number of ions is remarkably reduced from 54 to 16 
because of the extremely high computational costs. Over 1000 electronic states has to be 
included to describe a warm dense KC1 system consisting of 54 ions at several eV in the 
DFT calculation. It is challenging to perform such a calculation with all electronic states 
explicitly described at a plane wave cutoff of 400 Ry. 

To directly compare with the experimentally measured transition energy, an correction 
of 92 eV is added to the calculated edge position to compensate the intrinsic inaccuracy of 
the DFT method in calculating E\ s - ft. 


III. PRINCIPAL HUGONIOT AND THE INFLUENCE OF THE APPROXIMATE 
DESCRIPTION OF IONIZATION 


Calculated principal Hugoniot of KC1 is presented in Fig. UK a) together with available 


experimental measurements[3QH32], Theoretical principal Hugoniots derived from t 


model 


re QEOS 

12] (calculated using the mpqeos program included in the MULTI package 41]) as well 
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as Altshuler’s model 


29] are also displayed for comparison. Corresponding temperature along 


the Hugoniot is displayed in Fig. QJb). Fig. [Q(a) shows that the FPMD Hugoniot generally 
agrees well with experimental data. A small overestimation about 25 GPa in pressure is 


observec 
et al. 30 


at densities higher than 5 g/cm 3 in comparison with the measurements of Kormer 


The origin of this deviation is not clear so far. On the theoretical side, it is possibly the 


reflection of the overestimation of ionization in the FPMD method. Calculation on LiF 


42] 


displayed a similar trend compared with Kormer et aids measurements. On the other hand, 
uncertainties in experimental measurements could also be a cause of the deviation. Fig. [l](a) 
shows that the deviation starts from a sharp discontinuity taking place around 5 g/cm 3 
in the measured Hugoniot. This discontinuity does not belong to any of the first-order 


transitions recognized so far, as have been summarized by Duvall and Graham in Ref. 


43 


The corresponding pressure (150 GPa) of the discontinuity is far beyond the pressure of the 
polymorphic B1-B2 phase transition (~2 GPa[32]) and the pressure of the solid-liquid phase 
transition ( 33 GPa ~ 48 GPa). 43] 

The QEOS Hugoniot lies above the experimental measurements as well as the FPMD 
Hugoniot. The overestimation is derived from a much higher estimation to the temperature, 
as displayed in Fig. [0(b). It is quite unusual that the QEOS has a poor performance on a 
typical ionic crystalline material in the calculation of Hugoniot. It works reasonably well for 
other typical ionic crystalline materials like NaCl[l3] and LiF 42]. The Altshuler’s model 
has a similar performance as the QEOS model in the calculation of pressure. However, it 
gives an even higher estimation to the temperature and tends to overestimate the pressure 
on further compression. The large variation of temperatures displayed in Fig. [0(b) highlights 
the necessity to take temperature into account as a critical reference in building theoretical 
EOS models. 

Calculated degree of ionization along the Hugoniot is displayed in Fig. [2(b). Using the 
FPMD method, the statistics of electron is described by Mermin’s functional 2d| , which 
assigns an average occupation number to each electronic state accordin g to the Fermi-Dirac 
distribution. By this method, the ionization ratio R{T ) is estimated as 44] 


r 


D(E)dE 

e (E-/j,)/kBT + ’ 


(3) 


where A* is the total number of ions, Z is the averaged charge number for each ion, D(E) 
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FIG. 2. (Color online) (a)Energy gap E g and (b) ionization ratio along the principal Hugoniot 
calculated using the FPMD simulations. 


is the the density of state, and /i is the chemical potential of the system. 


As a result of the mass-action law, 45] the ionization ratio R{T ) is roughly proportional 


to exp\-(E g /2k B T)] at modest temperatures, where E g is the energy of the gap. For tem¬ 
perature much higher than Eg/2, /j would be much lower than the energy of electrons at the 
bottom of the shell, as the consequence of finite number of electrons. Accordingly, R(T ), 
as estimated by Eq. (J3J), would approach to 100%. For commonly used LDA and GGA 
exchange-correlation functionals in the FPMD simulations, E g is generally underestimated 
by 40% - 60% for nonmetals 15], which implies a sizable overestimation to the R(T) of KC1. 

The variation of E g along the Hugoniot is also an important factor affecting the estimation 
of R(T ). When temperature and density keep increasing along the Hugoniot, E g approaches 
to zero as the result of pressure ionization and thermal smearing. The vanishing of the band 
gap along the Hugoniot is broadly observed in materials of finite energy gap under ambient 


conditions, among which LiF 


42] is a recent example under extensive studies. Fi 


displays the variation of band gap along the Hugoniot, determined through Eq. 


44], It 


shows that the energy gap vanishes around 6.3 g/cm 3 , where the 3 p electrons of Cl are 
hybridized with the 3d electrons of K. This value can also be derived from the variation of 
density of states (DOS), as displayed in Fig. [3](d). It provides a threshold indicating when 
the system becomes completely metallic and the inaccuracy caused by the gap becomes 
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FIG. 3. (Color online) K-XAS of Cl and corresponding PDOS for KC1. (a) XAS calculated under 
ambient conditions; and (b) XAS under typical warm dense conditions, (c) Corresponding PDOS 
for (a) and (b), where an energy shift of 2820 eV is applied to move the chemical potential to the 
vicinity of the origin point, (d) PDOS before and after the vanishing of energy gap . 


negligibly small. Beyond the threshold, the reliability of the FPMD method has been well 
demonstrated by the calculations on typical metallic materials, such as A1 and Fe 


20 


22 , 


34, [46]. An accuracy of several percentages in both pressure and electronic structures can 


be achieved in the warm dense region (before the localized electrons in the lower shell are 
thermally excited). These results imply that the inaccuracy of ionization in the FPMD 
method is constrained by the pressure ionization and thermal smearing, and it will not 


accumulate globally, i.e., increase unboundedly, along t 


le Hugoniot. This is, in indeed, what 


displays in Fig. Ufa) and in the previous calculations 16-22], Extending this observation to a 
much higher temperature region (above 100 eV) strongly depends on the temperature effec 


of exchange-correlation functionals. Recent progress on this topic can be found elsewhere 


13 


47j. 


To quantitatively estimate how an approximate description of ionization in FPMD affects 
the Hugoniot result, we manually adjust the temperature of electrons. Besides the Hugoniot 
calculated at T e = T), Fig. [Ha) also displays the Hugoniots of T e = 300 K, and T e = (5/9)T/ 
The fraction of 5/9 gives an approximate correction to the overestimation of ionization 
caused by the underestimated band gap for the p < 4.0 g/crn 3 (T e < E g /2) part in the 
Hugoniot curve. The calculation at T e = T t gives a better estimation to the pressure for the 
p > 6.3 g/cm 3 part, when the energy gap vanishes and the material is completely metallic. 
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TABLE I. Contributions of ionic kinetic pressure electronic kinetic pressure P ey ki and virial 
pressure P po t to the total pressure (in GPa) at selected electronic temperature and ionic temperature 
for thermal states p = 3.9 g/cm 3 and 4.9 g/cm 3 corresponding to typical ionization states along 
the principal Hugoniot. 




3.9 g/cm 3 



4.9 g/cm 3 



P i,k 

P e,k 

P pot 

total P 

P e,k 

Ppot total 

H 

II 

EH 

5.07 

1262.82 

-1204.83 

63.06 



Tj=0.5eV T e =(5/9)Tj 

5.04 

1268.42 

-1205.98 

67.47 



T e =300K 

5.06 

1267.07 

-1205.67 

66.46 



Te=T i 

10.10 

1281.90 

-1210.36 

81.64 



Tj=1.0eV T e =(5/9)Tj 

10.06 

1279.66 

-1209.00 

80.72 



T e =300K 

10.10 

1280.92 

-1208.90 

82.12 



Te=T i 

20.19 

1305.98 

-1218.50 107.67 25.37 

1699.98 

-1547.22 178.13 

Tj=2.0eV T e =(5/9)Tj 

20.20 

1299.31 

-1214.45 

105.06 25.51 

1694.25 

-1541.46 178.29 

T e =300K 

20.22 

1296.38 

-1212.94 103.66 25.39 

1691.45 

-1538.74 178.10 

II 

<P 

H 




38.02 

1735.11 

-1563.49 209.63 

Tj=3.0eV T e =(5/9)Tj 




38.09 

1714.11 

-1550.48 201.72 

T e =300K 




37.98 

1705.38 

-1543.42 199.94 

II 

<P 

EH 




50.81 

1779.85 

-1582.51 248.16 

Tj=4.0eV T e =(5/9)Tj 




50.78 

1735.46 

-1558.39 227.86 

T e =300K 




50.72 

1722.14 

-1547.44 225.42 


Using these two calculations as references, the deviation of the calculation can be estimated 
as the smallest difference from the T e = Ti or T e = (5/9 )Ti curves to the experimental 
measurements. As displayed in Fig. [TJ the largest deviation in the calculation takes place in 
the region between p = 4 g/cm 3 and p = 6 g/cm 3 . Below that, the T e = 5/9 Ti curve only 
has a small correction of 1 ~ 2 GPa to the T e = T % curve because of the small magnitude 
of R(T), as displayed in Fig. [21(b). The largest deviation is about 50 GPa (~ 25 %) at 
p ~ 5.3g/cm 3 , including experimental uncertainties. This gives an estimation to the worst 
scenario in the Hugoniot calculation. 

To clarify the underlying mechanism of the small influence of the approximate ionization, 
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a detailed decomposition of pressure is given in Table U for selected electronic and ionic 
temperatures, corresponding to Hugoniot points at p = 3.9 and 4.9g/cm 3 . As displayed 
in Fig. [2^a), they represent the starting and ending points of the sharp transition of the 
band gap along the Hugoniot. In Table |H the pressure is divided into three parts according 
their origins as P(T, p) = P e ^(T, p ) + P i k {T, p) + P pot {T , p). Subscripts e, k and i, k represent 
kinetic pressures of electrons and ions respectively, and subscript pot refers to virial pressures 
contributed by ion-ion, electron-electron, and ion-electron interactions. According to the 
virial theorem 15], P k is proportional to 2 E k , and P po t is approximately proportional to 


Ppot- 

Table [I] displays the change in P pot and P £)k brought about by different ionizations. Gen¬ 
erally, P £) k grows with the increase of ionization, whereas P pot goes to the opposite direction. 
The fluctuation at low T e ( T e < 0.6 eV in the table) is caused by electronic excitations 
from localized electronic states to non-localized states near chemical potential. This is a 
characteristic of ionic materials. For KC1, the fluctuation is caused by the excitation of the 
localized 3 p electrons of Cl to the non-localized 4s and 3d electrons of K. Most of the changes 
are brought by the variation of P e k . The direction of the change depends on the competition 
between the average kinetic energies of the localized and non-localized electrons. 

The cancellation between the increase of P etk and decrease of P pot leads to a much smaller 
change in the total pressure. For the Hugoniot point at p = 3.9 g/cm 3 ( T t ~ 1.0 eV), the 
T e = (5/9)Tj calculation gives a better estimation to the ionization. Its correction to the 
total pressure is less than 1 GPa, which is about 1% of the FPMD result. For the Hugoniot 
point at p = 4.9 g/ern 3 ( 7) ~ 2.0 eV) the correction calculated with T e = (5/9)7) is less than 
0.5 GPa. The real correction should be much smaller than that, because the energy gap is 
now close to zero, as displayed in Fig. [2] (a). 


IV. K-SHELL XAS OF Cl AND EVOLUTION OF ELECTRONIC STRUCTURES 

XAS provide information on both ionization and electronic structures. Owing to the 
thermal excitation of electrons in WDM, electronic structures 10-20 eV below the chemical 
potential can also be detected by XAS. This is quite different from the low temperature cases, 
where only electronic structures above or close to the chemical potential is available for X- 
ray transitions. Accordingly, terminologies commonly used in XAS of condensed matter 
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FIG. 4. (Color online) Probability distribution of the 3 p wave functions of KC1 at E = 2816 eV. 
(a) The 3 p wave functions of Cl under ambient conditions, (b) The hybrid 3 p wave functions of 
Cl and K at 5.8 g/cm 3 . The K and Cl ions are represented by purple and green balls, respectively. 
The isovalue is set to be 0.001. 


should be used cautiously to avoid discrepancies. Here, we do not distinguish sub-categories 
of XAS and focus on absorption properties near the beginning of K-absorption spectra ( 
usually called the edge region) for Cl ions in the shock-generated warm dense KC1. Note 
that, the absorption edge in WDM could be far below (or above in some extreme cases) the 
chemical potential, as will be illustrated in the following part of this section. 

The K-absorption spectra under ambient conditions are compared with experimental 


measurements 


48] in Fig.[3](a). It serves as a reference illustrating how the calculated spectra 


can be interpreted. Fig. |3](a) shows that the calculation can reproduce the general feature 
of the spectra with small difference in details. The energy gap issue is paid special attention 
in the calculation. A correction of 3.7 eV obtained from a separate GoWo calculation[49] 
is added to the FPMD energy gap (~ 5 eV), so that the calculated gap value agrees with 
the experimental measurement of 8.69 eV [ll|. Hybrid functional HSE06 50] only slightly 
improves the gap to about 6 eV, and therefore is not helpful to the gap issue. Note that the 
GoWo correction is not included in the XAS calculation of warm dense states, because the 
accuracy of the method relies on the description of ionization in the FPMD method, which 
is generally overestimated before the vanishing of the gap. Some kind of self-consistent GW 
approach 51] is necessary to get a theoretically consistent result. 

Fig. [3](b) displays calculated K-absorption spectra of warm dense KC1. The prominent 
feature of the spectra is the persistence of the energy gap (around 2820 eV) under compres¬ 
sion, disappearing at a high density ~ 6.3 g/cm 3 . The figure also shows that the position 
of the chemical potential, around 2820 eV, is insensitive to the variation of density and 
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FIG. 5. (Color online) Transmission spectra of warm dense KC1 (a) at fixed density p = 8.3 g/cm 3 ; 
and (b) at fixed temperature T = 3.0 eV. 


temperature, which is similar to what was found in warm dense A1 34]. 

Projected density of states (PDOS) in Fig. |3](c) and (d) show that the electronic states 
across the gap have different origins. The majority of the electronic states below the gap 
consist of the 3 p states of Cl and K; whereas those above the gap mainly come from the 4s 
and 3d states of K. The figures also show a strong hybridization between the 3 p orbitals of 
Cl and K under compression. This trend is further illustrated in Fig. |3](b) by the increase 
of band width in the sub-band below the gap (below 2820 eV). According to the theory 
of pressure ionization 52, [53], this hybridization of electronic states of different ions is an 
evidence of the pressure ionization, i.e., a transition from localized states to non-localized 
states. This transition is also confirmed by the variation of wave functions of the 3 p orbitals 
at E — 2816 eV, as displayed in Fig. [4] for warm dense states p = 2.0 g/cm 3 (T = 300K) 
and p = 5.8 g/cm 3 (T = 3.3 eV). Here, 2816 eV is the energy of the maximum PDOS for 
the 3 p states of Cl under ambient conditions, and the thermal state of p = 5.8 g/cm 3 is the 
last state in the calculation just before the closing of the gap. The PDOS after the closing 
of gap is also displayed in Fig. [3](d). These results show that the 3 p states of both Cl and 
K are part of the continuum (non-localized states), and the gap around E = 2820 eV is in 
the middle of unbounded states above the bottom of the continuum. 

Besides the absorption cross section a, the transmission (or attenuation) of the material 
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Q 


was also measured in experiments. On many occasions, it was used to determine the ab¬ 
sorption edge positions)^, 101. The transmission is proportional to exp(-crh), where h is the 
areal number density. Fig. 0 displays the dependency of transmission spectra (corresponding 
to the K-absorption of Cl) on density and temperature. The h here takes the value of 790 
A -2 , derived from recent experiment of Zhao it al. 8] Fig. 0(a) displays the variation of 
transmission as a function of temperature, at a fixed density of 8.3 g/cm 3 . It shows that the 
absorption edge decreases with the growth of temperature, which is mainly caused by the 
thermal broadening of depopulated electronic states below the chemical potential. Fig. 0(b) 
presents the transmission spectra at fixed temperature of 3 eV, showing that the absorption 
edge is insensitive to the variation of density. Similar trend has also been found in warm 
dense Al, 34] indicating that the K-absorption edge is useful to diagnose the temperature of 
warm dense state for a variety of materials. The feature corresponding to the energy gap 
appears as a peak around 2820 eV in Fig. 0 which gives an explanation to the “bump” 


feature observed in the time-resolved transmission spectra measured by Bradley et al. 


io| 


To compare with the experiments of Zhao et al, jg] in which the intensity of transmitted 
X-ray was measured, the intensity 1(E) is further derived from a according to 1(E) = 
Io(E) exp[-a(E)h], with Iq(E) the intensity of backlighter. Fig. 0 supplies a comparison 
between calculated and measured 1(E). The background is deducted according to the left 
part of the measured spectra below the absorption edge. The intensity of backlighter jg] 
measured in a separate experiment is also included as references, and the energy scale is 
shifted by 2820 eV so that the position of the chemical potential is located at the vicinity 
of the origin point. Fig. 0(a) displays the calculation of KC1 under ambient conditions. 
It well reproduces the XAS up to ~30 eV above the absorption edge. Fig. 0(b) and (c) 
show that our calculation is in reasonable agreement with the experimental measurements 
up to the chemical potential. For the next 20 eV above that, the calculated intensity is 
lower than the experimental measurements by about 50%, which suggests that a refined 
theoretical treatment is necessary to deal with the spectra above the energy gap in warm 
dense materials. 
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FIG. 6. (Color online) Transmitted X-ray Intensity of KC1 for (a) cold samples; (b) 10.5 g/cm 3 ; 
and (c) 8.3 g/cm 3 . Experimental measurements and intensity of backlight are also presented for 
comparison. The energy scale is shifted by 2820 eV to put the chemical potential at the vicinity 
of the origin point. 


V. SUMMARY 


In summary, we investigate the variation of electronic structures and the influence of the 
approximate description of ionization to the calculation of pressure through first-principles 
calculation of principal Hugoniot and XAS for warm dense KC1. We show that pressure ion¬ 
ization and thermal smearing are the major factors to limit the deviation in the calculation 
of pressure. They prevent global accumulation of such deviations along the Hugoniot. In 
addition, the cancellation between kinetic pressure of electrons and virial pressure of inter¬ 
actions further reduces the error. The calculation of XAS shows that the band gap of KC1 
persists after the occurrence of the pressure ionization of the 3 p electrons of Cl and K, which 
provide a detailed understanding to the evolution of electronic structures for WDM. 
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